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Abstract 

The next generation of long baseUne neutrino experiments will aim at determining the value 
of the unknown mixing angle, ^13, the type of neutrino mass hierarchy and the presence of CP- 
violation in the lepton sector. Beta-beams and electron capture experiments have been studied 
as viable candidates for long baseline experiments. They use a very clean electron neutrino beam 
from the /3-decays or electron capture decays of boosted ions. In the present article we consider 
an hybrid setup which combines a beta-beam with an electron capture beam by using boosted 
Ytterbium ions. We study the sensitivity to the CP-violating phase 6 and the ^13 angle, the CP- 
discovery potential and the reach to determine the type of neutrino mass hierarchy for this type 
of long baseline experiment. The analysis is performed for different neutrino beam energies and 
baselines. Finally, we also discuss how the results would change if a better knowledge of some of 
the assumed parameters was achieved by the time this experiment could take place. 
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I. INTRODUCTION 

In the past decade, atmospheric solar 0,Q,fl, reactor and long-basehne 

accelerator j^, neutrino experiments have provided compelling evidence for the phe- 
nomenon of neutrino oscillations. This has reshaped our understanding of the properties of 
elementary particles as it implies that neutrinos have mass and mix. The combined data 
can be described by two mass squared differences, Amfi and Am^i, where Am^j = m'j — , 
whose current best fit values are |Am|i| = 2.4 X 10-3 eV^ and Am^^ = 7.65 x 10"^ eV^ 
The two mixing angles 612 and 623 drive the solar and KamLAND, and atmospheric and 
MINOS neutrino oscillations, respectively, and are measured to be sin^ 612 = 0.304 and 



sin^ 623 = 0.50 



The third mixing angle, 613, is yet undetermined but is known to be 



small or zero. With available data, 613 is constrained to be 



sin^ ^13 < 0.040 (0.056) at 2a {3a) . (1.1) 



It is interesting to note that very recently a first hint in favour of ^13 7^ has been found 12 



in a combined analysis of atmospheric, solar and long-baseline reactor neutrino data, with: 

sin^ 013 = 0.016 ±0.010 at la , (1.2) 



implying a preference for 613 > at 90%CL. A different analysis [13|] confirms the hint for 
6*13 7^ at 1.5 cr from the analysis of solar and KamLAND data owing to the latest SNO 
results, but not the one from the atmospheric data. 

Although the experimental progress in neutrino physics over the last decade has been 
conspicuous, many of the fundamental questions surrounding neutrinos still need to be 
addressed. Understanding of the physics beyond the Standard Model responsible for neutrino 
masses and mixing requires knowledge of the nature of neutrinos (whether Dirac or Majorana 
particles), the neutrino mass ordering (normal or inverted), the absolute neutrino mass scale, 
the value of the unknown mixing angle ^13, and whether CP-symmetry is violated in the 
lepton sector. It will also be necessary to improve the precision on the known parameters, 
in particular to measure any deviation from maximal 623 and, if so, to determine its octant. 

Some of the issues above will be addressed by a future program of neutrino oscillation 
experiments 1^, [l^ . In particular, long baseline experiments using conventional beams [l^ 
and nuclear reactors |l^ will be the first to explore 613 below the current limit and maybe 
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confirm the hint for ^13 7^ [12]. If ^13 is close to the present bound imposed by the running 



and near future experiments, the next generation of superbeams |17l. |18[|. an extension of a 
conventional beam with an upgrade in intensity and detector size, and wide-band beams 
will probe CP-violation and, for sufficiently long baseline, the neutrino mass hierarchy. For 
small values of 61^ or, if ^13 is large but a better precision on the neutrino parameters needs to 
be achieved, the community must turn to the novel concepts of the neutrino factory 20|, l2l| or 



beta-beam 



22 



23| . Whereas conventional beams sourced from pion decays have an intrinsic 



contamination of electron neutrino at the ~ 1% level (owing to kaons in the beam), neutrino 
factories and beta-beams will have clean sources from highly accelerated muons and ions, 
respectively, producing a well-collimated beam. In a neutrino factory, muons (antimuons) 
are produced, cooled and accelerated to a high boost before being stored in a decay ring. The 
subsequent decay sources a muon neutrino (muon antineutrino) and electron antineutrino 
(electron neutrino) which are aimed at magnetised detectors located a very long distance 
from the source. The use of magnetised detectors is necessary to separate the 'right muon' 
disappearance signal from the 'wrong muon' appearance signal, which is sensitive to matter 
effects and CP-violation. A beta-beam will exploit accelerated ions that /?-decay sourcing 
a clean, collimated, electron neutrino beam. Magnetised detectors will not be necessary in 
this case, the only requirement being possession of good muon identification to detect the 
appearance channels. Therefore, water Cerenkov (WC), totally active scintillator, liquid 
argon detectors and non-magnetised iron calorimeters could be used, depending on the peak 
energy. 

The determination of the oscillation parameters is severely affected by degenera- 

the possibility that different sets of the unknown parameters 



cies m 
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26 



271, 



(sgn(Am|]^), 5, ^13, $23 octant) can provide an equally good fit to the probability for neutrino 
and antineutrino oscillations, for fixed baselines and energy. Therefore, a high precision mea- 
surement of the appearance probabilities is not sufficient to discriminate the various allowed 



solutions. In order to weaken or resolve this issue, various strategies have 



)een put forward: 



exploiting the energy dependence of the signal in the same experiment [19|, |29| , using reactor 
neutrino experiments with an intermediate baseline 3^, combining different long baseline 
experiments [si], adding the information on ^lofrom reactor experiments [32], or using more 



than one baseline for the same beam 



33, 



34, 



Earth matter effects in multi-GeV atmospheric 39 



36 



37 



40 



41 



38 
421 



In addition, 613 controls the 



431] and in supernova neutrino 
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oscillations [443 (^^^ ^^^^ l^)- These might provide useful information on the type of 
neutrino mass hierarchy and 6*13; the magnitude of the T- violating and CP- violating terms 
in neutrino oscillation probabilities is directly proportional to sin ^13 



46 



47|. 



In beta-beam experiments, the energy dependence of the signal is typically used to extract 



information on the mass hierarchy and CP- violation. Matter effects increase with base 
and energy suggesting that setups with baselines > 600 km are necessary 



29 



38 
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me 



50 



51 



49, 



52| for the determination of the type of neutrino mass ordering. Such strategies 



would make use of a proposed upgrade to the CERN Super Proton Synchrotron (SPS) 
which would equip the accelerator with fast superconducting magnets allowing high boosts 
and fast ramps. The latter are important to reduce the loss of ions through decay in the 
acceleration stage. A sister approach to the beta-beam is to use the neutrinos sourced from 
ions that decay mainly through electron capture |53|, |5J, l55|, l56|- If the electron capture 



decay is dominated by a single channel, then a monoenergetic electron neutrino beam can be 
produced this way. In this case, all the beam intensity can be concentrated at the appropriate 
energy to get the best sensitivity to the oscillation parameters. In order to disentangle the CP 
violating phase with neutrino beams only, one makes use of the different energy dependence 
of the CP-even and CP-odd terms in the appearance probability [47]. Electron capture 
competes with /S^-decay when the Qec- value > 2me, rrie being the electron mass. With the 



ions identified in 



531 ] ■ the use of an upgraded SPS or the Tevatron ^ is necessary to source 



baselines in excess of CERN-Frejus (130 km). In this paper, we discuss a hybrid of these 
two approaches. ^ By selecting a nuclide with Qec ~ 4 MeV, we can make use of neutrinos 
from an electron capture 5-spike and continuous spectrum simultaneously. Assuming a 
detector with low energy threshold, the use of such ions allows one to exploit the information 
from the first and second oscillation maxima with a single beam, in a similar way to the 
approach used in Ref. |29|. There the spectral information was used to remove some of 
the degeneracies and reach physics sensitivities comparable to the scenarios with a neutrino 
and antineutrino beam that are often presented in the literature. The use of the hybrid 



^ Note that the present Tevatron configuration does not ramp fast enough, resulting in a high loss of ions, 

so this might not be a very realistic experimental setup, at least in thejjresent configuration. 
^ Note that the use of this hybrid approach was first mentioned in Ref. 55| , although the proposal, unlike 
the present case, was to use long-lived ions. In addition, the phenomenology of this approach was not 
studied. 
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approach we propose makes it possible to use a monochromatic beam at higher energies and 
a beta-beam at lower energies. The need for good neutrino energy resolution at the higher 
energies will therefore be less crucial than for high-7 beta-beam scenarios. 

The paper is organised as follows. In Sec. [TTl we introduce the hybrid idea and show 
two possible ions which could be used. In Sec. [1111 we present the different set-ups studied, 
discussing the choice of boost factor, basehne and various detector options. The description 
of the simulation and analysis is presented in Sec. IIVI and the results in Sec. |Vl The final 
discussion and conclusions are drawn in Sees. |Vl] and I VI II 



II. THE BETA-BEAM AND ELECTRON CAPTURE COMBINATION 



In this section we introduce the idea of the beta-beam and electron capture hybrid ap- 
proach. We present the spectra of the two branches, their ratio and discuss two nuclides 
which have desirable properties. 



The beta-beam is a proposal, originally put forward by P. Zucchelli 22|, to accelerate 
and then store /3-emitting ions, which subsequently decay to produce a well coUimated, 
uncontaminated, electron neutrino (or antineutrino) beam. The high luminosities required 
to achieve a useful physics reach point towards ions with small proton numbers to minimise 
space charge and half-lives ~ 1 second to reduce ion losses during the acceleration stage 
whilst maintaining a large number of useful decays per year. The most promising candidate 
ions are ^^Ne and for neutrinos, and ^He and ^Li for antineutrinos. A variant on the 
beta-beam idea is the use of electron capture to produce monoenergetic neutrino beams. 
Electron capture is the process in which an atomic electron is captured by a bound proton 
of the ion A{Z, N) leading to a nuclear state of the same atomic number A, but with the 
exchange of the proton by a neutron and the emission of an electron neutrino, 

A{Z,N) + e- ^ A{Z -1,N +1) + . (2.1) 

The idea of using this process in neutrino experiments was independently discussed in 



Refs. 



53, 



55l |. In Ref. |56l], ions with low Qec- value and long half-life, such as ^^''Sn, 
were proposed to be accelerated to very high boosts with the LHC. Baselines of 250 km 
and 600 km were considered with the spectral information coming from the position of the 
events in the detector. Sensitivities comparable to a Neutrino Factory were obtained for a 
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single boost. However, in order for electron capture machines to become operational, nuclei 
with shorter half-life are required. The recent discovery of nuclei far from the stability line 
with kinematically accessible super-allowed spin-isospin transitions to giant Gamow- Teller 
resonances (see, for example, Ref. [s^]) opens up such a possibility. The rare-Earth nuclei 
above ^^^Gd have a short enough half-life to allow electron capture processes in the decay 



ring, in contrast to fully-stripped long-lived ions 55|, l56| • This was the alternative put for- 



ward in Ref. 53| where the use of short-lived ions with QEC-values around 1-4 MeV was 
proposed. Machines such as the SPS, an upgraded SPS and the Tevatron could then be 
used for the acceleration. The ion ^^°Dy, with Qec- value 1.4 MeV, was investigated for the 
CERN-Frejus (130 km) and CERN-Canfranc (650 km) baselines and different boost factors. 
It was found to have very good physics reach 0, [s^ . Owing to the monochromatic nature 
of the beam, multiple boosts are necessary to resolve the intrinsic degeneracy in this case. 

In the following, we demonstrate how the flux for the electron capture/beta-beam can 
be built up by discussing them separately and comparing branching ratios. Let the mass 
difference between the parent and the daughter nuclei, AM^ = Ma{Z,N) — Ma{Z — 
1, + 1), include the mass and the binding energy of an atomic electron as well. For 
electron capture, the maximum kinetic energy release is thus given by Qec = AM^ . For 
/3"'"-decay, however, the final nucleus has an excess electron since a positron is produced. 
The maximum kinetic energy release is thus given by Qp+ = AM^ - 2 m^. Clearly for 
(AM^ =) Qec < 2me, electron capture is the only allowed process for a proton-rich nucleus. 
For Qec > 2me, electron capture and positron emission compete, their branching ratios 
dependent on Qec- If decay through a emission is also allowed, it is important that this has 
a relatively low Q- value so as not to be the dominant channel ^. For a number of useful ion 



decays per year Anions, the electron capture neutrino flux is given by 



53 



54| 



for each decay channel. Here, L is the baseline, 7 is the Lorentz boost, i?^^ (= Qec) is the 
neutrino energy in the ion rest frame and Ei, is the neutrino energy in the lab frame. 

The flux for the /5-spectrum is found in the usual way. In the rest frame of the ion, the 



^ The a decay branching ratio is strongly dependent on the QEC-value. For low QeCj the a decay probability 
is sufficiently long as to allow the weak decay modes to be the main channels. 
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electron neutrino flux is proportional to 



(2.3) 



Here, E. 



d cos OdEj^i 

+ me = Qec ~''^e) is the total end-point energy of the decay. The neutrino 



flux per solid angle at the detector located at distance L from the source after boost 7 is 



48| 



dVtdy 



Anions 7' 



(2.4) 



where < y 



2^ 



<l-ye,ye = m^/E^, and 



60 




y^,{2-9yi-8yt) + 15yt\og 



1 



1 



Ve 



(2.5) 



Similarly to the case of electron capture, a neutrino with energy E^.^ in the rest frame will 
have a corresponding energy E,y = 2'-)E^i in the laboratory frame along the ^ = 0° axis. 

All the known nuclear structure information on the A = 148 and A = 156 nuclides has 
been reviewed in Ref. 



and Ref. 



591], respectively, where the information obtained in 



various reaction and decay experiments is presented, together with adopted level schemes. 
Currently, a systematic study of electron capture decays^ in the region of ^^^Gd, relevant for 
monoenergetic neutrino beams, is being carried out [60]. Here, we consider two nuclides, 
7Q^Yb and gf^'^Tb, that decay through electron capture and /?"'"-decay with similar branching 
ratios whose lifetimes are not too long or too short. Their decays are summarised in Tables [T] 
and [TTl Ytterbium is a nuclide 70^ Yb with spin-parity 0^, which decays 90% via electron 
capture plus (3~^-decay l^^, with 38% via electron capture and 52% via (3~^-decay 60] . The 
remaining 10% goes into alpha particles and a different flnal state. This relatively small 
branching ratio into alphas helps the nuclide to have a short enough half-life, 26.1 seconds. 
It is important to note that this electron capture-/5"'"-decay transition has only one possible 
daughter state with spin-parity l"*", i.e., it is a Gamow- Teller transition into an excited state 
of Thulium, gg^Tm*. The transition QEC-value is Qec- value = 3.58 MeV. However, the 
excitation energy of the flnal nuclear state (0.12 MeV) needs to be taken into account and 
thus, the effective QEC-value (difference in the total kinetic energies of the system after and 
before the decay) is 3.46 MeV jsoj. The electron capture energy of ~ 4 MeV is well suited 



QEC-values are typically calculated between ground states unless stated otherwise. 
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Decay 


Daughter 


Neutrino energy (MeV) 


BR 


/?+ 
EC 
a 


68 ^1 


2.44 (endpoint) 

3.46 

4.81 


52% 
38% 
10% 



TABLE I: Decay summary for 7o^Yb. The QEC-value for the transition between ground states is 
3.58 MeV and taking into account the excitation energy of the final nuclear state (0.12 MeV), the 



effective Q|f;- value is 3.46 MeV 



59 



Decay 


Daughter 


Neutrino energy (MeV) 


BR 


EC 


148p J* 
64 

^fCd* 


2.05 (endpoint) 
3.07 


32% 
68% 



TABLE IL Decay summary for gg^^Tb. The QEC-value for the transition between ground states 



is 5.77 MeV and the effective Qgc-value to the excited state is 3.07 MeV 



58|, 



61 



62| 



to the intermediate-baselines of Europe and the USA with the available technology, or those 
available with future upgrades. On the other hand, the gf^'^Tb isomer with spin-parity O"*" 
has a Qsc-value of 5.77 MeV 



58l . l6ll ]. Although the decay to the ground state of gl^Gd is 
highly forbidden, the presence of a Gamow- Teller resonance allows the decay into an excited 
state with effective Q- value 3.07 MeV 62|. This nuclide is longer lived than yg^Yb (with 
a half-life of 2.2 minutes) and will require slightly higher boosts. It is still well suited to 
intermediate baselines. However, the dominance of the electron capture over the /5^-decay 
channel makes this nuclide less desirable. The count rate will be dominated by the single 
energy of the electron capture which provides insufficient information to obtain t he g ood 
sensitivities aspired to by future long baseline experiments. It was shown in Refs. [ssl. Is^ 
that two runs with different boosts are necessary for an exclusive or dominant electron 
capture channel to break the intrinsic degeneracy and achieve good CP-violation discovery. 
Hence, in what follows we will study this hybrid approach focusing on ^^^Yb. 
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III. EXPERIMENTAL SETUPS 



In this section we consider different boosts, baselines and detectors. We first discuss the 
available (or possible future) accelerator technology and identify the possible boost factors 
in combination with different baselines. We then discuss the main characteristics of the 
detectors considered in the analysis. 



A. Choice of 7 and baseline 

In this paper we consider the use of a neutrino beam sourced from boosted ^^^Yb ions 
directed along a single basehne. As described above, both the electron capture and /3+-decay 
channels are to an excited state of ^^^Tm with a QEC-value of 3.46 MeV. In order to fully 
exploit the electron capture deca^ mode, the nuclides cannot be fully stripped; at least 16 
electrons being left on the ion |[631]. The maximum boost, 7max, available is thus 

Eacc Z -16 

7ma^ = , (3.1) 

mp A 

where rrip is the mass of the proton and Eacc is the maximum energy accessible with the 
accelerator. Current and future accelerator facilities would be an ideal production environ- 
ment. In this analysis, we consider the maximum boosts available from the current SPS and 
upgraded SPS (see Table IIIII) for the following baselines: 

1. Boost 7 = 166 with current SPS 

• CERN-Frejus (130 km) 

• CERN-Canfranc (650 km) 

2. Boost 7 = 369 with an upgraded SPS 

• CERN-Canfranc (650 km) 

• CERN-Boulby (1050 km) 

With the current magnetic rigidity of the SPS, the electron capture spike can be placed 
on first oscillation for the CERN-Canfranc baseline (650 km) with the beta-beam spectrum 
peaking around the second oscillation maximum (see Fig. [T]). 
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0.04 




(GeV) 




FIG. 1: Top panel: appearance probabilities for the CERN-Frejus (130 km) and CERN- 

Canfranc (650 km) baselines. The unoscillated t'e flux in the laboratory frame is shown for ^^^Yb 
given a boost 7 = 166 in arbitrary units. Bottom panel: appearance probabilities for the 

CERN-Canfranc (650 km) and CERN-Boulby (1050 km) baselines. The flux from a boost 7 = 369 
is shown in arbitrary units. In both cases, the blue lines correspond to CERN-Canfranc; the red 
being CERN-Frejus (top panel) and CERN-Boubly (bottom panel). The solid lines correspond to 
5 = 0°, dashed 6 = 90° and dotted 5 = —90°. The value sin^ 2^13 = 0.01 was taken for all curves. 
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Machine 


Tmax 


27maxQlfc (GeV) 


27maxQ^'l (GeV) 


SPS 


166 


1.15 


0.81 


Upgraded SPS 


369 


2.55 


1.80 



TABLE III: Maximum boosts and neutrino endpoint energies for ^^^Yb available for the current 
SPS setup and a proposed 1 GeV upgraded SPS. 

A detector with a low energy threshold is necessary to exploit the oscillatory structure of 
the appearance probability. The second option would make use of the upgrades to the CERN 
accelerator facilities necessary for suggested LHC upgrades. With a 1 TeV SPS, the electron 
capture beam could be placed at first oscillation maximum, or on the probability tail, for the 
CERN-Boulby baseline (1050 km) or, using the Tevatron, for the FNAL-Homestake baseline 
(1280 km). With an energy threshold of 250 MeV, these setups could exploit the information 
at second oscillation maximum to resolve some of the degeneracies in an approach analogous 
to Ref. 291, where the boost for ^'^Ne was chosen so that the boosted spectrum covered both 
first and second maximum. Owing to the spectral nature of the decay, in the analysis in 
Ref. 29| it was difficult to determine to what extent the highest energies contribute to the 
overall sensitivity of the setup. With the electron capture beta-beam combination we can 
investigate this issue as it is possible to place the electron capture spike on first oscillation 
maximum or on its tail whilst the beta-beam spectrum has minimal coverage of that energy 
range. In our analysis we perform {613,6) sensitivity contours through the consideration 
of each decay channel separately and their combination. This allows us to investigate the 
importance of the contribution of each decay channel and evaluate if the importance of the 
low and high energy contributions to the overall sensitivity. 



B. Choice of detector 



In a beta-beam, one aims to exploit the Ue z/^ channel; a detector with excellent 
muon identification capabilities and efficient neutral current background rejection is therefore 
required. For energies below ~1 GeV, water-Cerenkov detectors are typically chosen with 
the muons identified through the use of quasi-elastic events (QE). Efficient reduction of 
neutral current events and (subdominant) pions is through the identification of the decay 
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process. For higher energies, the number of QE events drops sharply where the deep inelastic 
scattering (DIS) component dominates the cross section. Water-Cerenkov detectors are 
usually not the best choice for the higher energies owing to high backgrounds and poor 
neutrino energy resolution. For large boosts, the 50 kton class of detectors such as liquid 
argon (LAr), using time projection chamber techniques, or total active scintillator detector 
(TASD), based on tracking calorimeter principles, are usually considered. In addition to QE 
events, these technologies also measure the energy deposited through the hadronic channels 
and DIS events are in principle also distinguishable. Their main disadvantage is their size 
which is far smaller than the next generation water- Cerenkov detectors discussed in the 
literature which typically have fiducial masses in the megaton scale. 

For a pure electron capture machine, the choice of detector technology does not depend 
on its energy reconstruction capabilities. In this case, the neutrino energy is given by the 
choice of ion and boost factor leaving no need to reconstruct the neutrino energy in the 
detector. For the hybrid approach we consider, it is possible, in principle to separate the 
energy of the line spectrum from the continuous spectrum. Suppose we identify an event 
and classify it as being a QE elastic event with energy E^^QE), then it must be the case that 
the true energy > E^i^QE). Thus, if one measures E^i^QE) > 2jEq, then this event 

must be attributed to the electron capture flux and hence, it is not necessary to reconstruct 
more precisely the true neutrino energy. The separation between the energy of the electron 
capture spike and the end-point energy of the beta-spectrum is 2me7. This should render 
the distinction between the electron capture and the beta-beam fluxes possible. We will 
assume this throughout the paper. 

In the analysis, we will follow two strategies regarding the detector type. We consider a 
generic detector technology, which could be LAr or TASD, with a flducial mass of 50 kton 
and assume that the neutrino spectral information can be extracted from the charged current 
events. On the other hand we also consider a 0.5 Mton (flducial) water-Cerenkov detector. 
In this case, following the prescription described above, we assume the neutrino energy 
from beta-beam events can only be reconstructed for QE events. However, we do include 
the information from the inelastic events. As no spectral information is possible for those 
events, we include them in a single bin. We will take perfect efficiency for the 50 kton 
detector (which can be easily scaled), and 70 % efficiency for the water-Cerenkov detector. 
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IV. SIMULATION AND ANALYSIS 



Based on the selected boost factors, baselines and type of detector we have discussed in 
detail in the previous section, in our analysis we will study and compare six different setups: 

1. 50 kton detector (LAr or TASD) with 2 x 10^^ ions/yr 

• Setup I: CERN-Frejus (130 km) and 7 = 166 

• Setup II: CERN-Canfranc (650 km) and 7 = 166 

• Setup III: CERN-Canfranc (650 km) and 7 = 369 

• Setup IV: CERN-Boulby (1050 km) and 7 = 369 

2. 0.5 Mton water-Cerenkov detector with 2 x 10^^ ions/yr 

• Setup III-WC: CERN-Canfranc (650 km) and 7 = 369 

• Setup IV- WC: CERN-Boulby (1050 km) and 7 = 369 

We will take a running time of 10 years for all the experimental configurations considered. 
The number of events is computed for each energy bin i, given by 



where A/" is a constant which takes into account the efficiency, the mass of the detector and 
the running time, ^{E) is the neutrino flux spectrum at the detector, P{E) is the probability 
function, cr{E) represents the total, QE or DIS cross section (as described in the previous 
section) and Ki{E) is the energy smearing kernel for the ith bin for which we take a Gaussian 
energy resolution function with a constant width of 150 MeV. 

In the experimental simulations performed in this study, our analysis is based on the 
following definition 



where are the predicted (or fitted) number of events for a certain oscillation hypothesis, 
and n\ are the simulated "data" for the true values of the parameters. The covariance matrix 
C given by 




(4.1) 




(4.2) 



Cij = 5y(5n-)^ , 



(4.3) 
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where {5n\) = y^n* + (/sys • njy, contains both statistical and a 2% overall systematic error 
(/sys = 0.02). In addition, we assume an intrinsic beam background of 0.1 % of the unoscil- 
lated spectrum originating from neutral current pion production and muons misidentified as 
electrons. In the energy range of interest, there are about 30 atmospheric neutrino events 
per kton-year which could mimic a muon coming from a ^ i'^ oscillation. We take 10~^ 
as the accelerator duty factor so that this atmospheric background amounts to 0.03 events 
per kton-year. 

In all setups considered in this paper, for the beta-beam part we take an energy threshold 
of 250 MeV and use 200 MeV bins above that value, except for non-QE events in water- 
Cerenkov detectors which are grouped in a single bin. For the electron capture events we 
always take a single bin. Unless otherwise stated, we impose restrictions on certain subsets 
of the fitted parameters in order to account for external information from other experiments. 
Thus, we introduce the so-called priors. Hence, if we want to restrict some parameter k, we 
introduce the central value Hc of the prior and the input error a^, and the actual minimisation 
is performed over the modified function 

X'--X'+ ^-^-^ . (4.4) 

(T 

In this work, we set priors for the experimentally known oscillation parameters, taking their 



best fit values as central values ll|] of the corresponding priors, and the half width of one 
standard deviation of their best fit values as the corresponding input errors. Specifically, 
Am2i, lAmgj^l, sin^ 612 and sin^ ^23 have been given errors 3%, 5%, 7% and 14% respectively. 
We marginalise over ^13 and 5 over their entire range. 

V. RESULTS 

In this section we present and discuss the results of our detailed numerical analysis of the 
various setups. In order to understand some of the features of these results, it is useful to 
consider an analytical approximation for the oscillation probability which for these energies 



and basehnes is given by [6J] (see also Ref. |65l|): 



^ - ' (A-Ai3)L 



P(z/e ^ z/^, L) ~ sin^ sin^ 2^13 (^^) sin^ 

+ cos ^13 sin 2^13 sin 2^23 sin 2^12 ^^rfl^ sin (f^) sin (^^^^^) cos + 5) (5.1) 
+ cos2^23sin2 2^i2(^)'sm2(^), 
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where A12 = Am2i/{2E) and A13 = Am'^i/{2E). We use the constant density approxima- 
tion for the index of refraction in matter A = y/2GFne{L), with ne{L) = 1/L j^^ n^{L')dL' 
the average electron number density. We analyse the sensitivity to ^13 and 5 for all the 
setups, and discuss the discovery reach for CP-violation and the type of neutrino mass 
ordering. 

A. Sensitivity to ^13 and 6 
Setups I and II 

Setups I and II use a low boost factor, 7 = 166 and a 50 kton detector, implying relatively 
low count rates. We find that the sensitivity to ^13 and 5 is very limited as a consequence. 
Violation of CP can be established only for a small range of values of the 5 phase and only 
if 6'i3 is close to the present bounds. The beta-beam channel contributes very little to the 
overall sensitivity of the setup. This is due to the 7^-dependence of a beta-beam flux. The 
small flux, when combined with the small cross-sections at the energies centred on second 
oscillation maximum, supplies a scarce count rate. The bulk of the sensitivity is therefore 
due to the electron capture channel placed on first oscillation maximum, as seen in Fig. [2l 
The performance of the beta-beam (upper row), electron capture (middle row) and their 
combination (lower row) is shown for 6'i3 = 1° and 5 = 90° for setup I (left column) and 
setup II (right column), at 90%, 95% and 99% confidence level (CL). 

Setups III and III-WC 

We have also examined the effect of placing the electron capture beam in the tail of first 
oscillation maximum. This gives the beta-beam coverage of the second oscillation maximum 
and substantial portions of the first oscillation maximum. For the CERN-Canfranc baseline 
and a boost 7 = 369, we find that the roles of electron capture and the beta-beam are 
reversed compared to 7 = 166. The beta-beam spectrum is peaked around ~ 1 GeV while 
the EC energy is of 2.55 GeV. As shown in Fig. [H the neutrino beta-beam spectrum is a good 
fit for the appearance probability at Canfranc, its peak sitting around the first oscillation 
maximum. The beta-beam now contributes much more to the sensitivity as it provides 
substantial information from the first oscillation maximum and a much higher count rate 
from the second oscillation maximum. 

In Fig. [3] and Fig. |4]we show the 90%, 95% and 99% CL contours for the setup III and 
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FIG. 2: 90%, 95% and 99% CL contours for setup I (left panels) and setup II (right panels). The 
parameters ^13 = 1° and 6 = 90° have been taken assuming normal mass ordering and 823 = 45°. 
The upper row is the contribution of the beta-beam, the middle row is that of the electron capture 
channel with the lower row being the combination. 
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setup III-WC, respectively, assuming that the hierarchy is normal and ^23 = 45°. In Fig. [3l 
we present contours for 5 = 90° and for the cases ^13 = 1° (left column) and 3° (right 
column). Similarly to Fig. [21 in Fig. [31 we have included the contributions from both the 
/3"'"-decay (upper row) and electron capture (middle row) channels separately to investigate 
their relative impact, in addition to the total sensitivity (lower row). 

Fig. [H shows the results for 61^ = 1° (left column) and 3° (right column) and four values 
for 6, including the hierarchy clone solution. Comparing the results for the setup III and 
setup III-WC, one can understand the effect of the event rate. We see that the 0.5 Mton WC 
detector gives a much better resolution, although not as much as one would naively expect 
from its larger size. As commented above, the substantially bigger size of the detector cannot 
be fully exploited due to the fast drop of the QE cross section at energies above 1 GeV, 
where the first oscillation maximum lies. 

The full power of the combination between the beta-beam spectrum and the EC channel 
is best illustrated in Fig. [3l We see that each of the two techniques separately suffer from a 
continuum of solutions. The shape of the allowed region in the 613 and 6 plane can be easily 
understood by looking at the form of the oscillation probability in Eq. (15. ip . Owing to the 
relatively short baseline and small matter effects, we can neglect for simplicity the matter 
potential term. A, in the following. As we measure the probability in only one polarity, for 
fixed energy and baseline, as is the case for the EC signal, we can relate the allowed values 
of 6 and 613 from a measurement of the probability as 



where is a constant which depends on the true values of 613, and 6. This approximation is 
valid as far as sin^ 29i3 ^ few x 10~^. The form of the expression in Eq. (15.21) matches the 
continuum of solutions in Fig. [3l In particular, we note the presence of a minimum for ^^13 at 
6 = — A13L/2. This result is more general then the approximated form of sin ^13 in Eq. (15. 2p 
and holds also for small values of sin 6*13. For the choice of the parameters used in Fig. [3l 
we have S ~ —44°. The range of the allowed solutions for 613 can be computed by looking 
at the amplitude of the function in Eq. 15.21 and is found to be A613 ~ ^ sin(A3iL/2) ' 
the choice of true values in Fig. [3l A6'i3 ~ 2°. 

Now consider the contribution of the beta-beam. For simplicity, we analyse its impact 
by looking at the energy of the first oscillation maximum, again neglecting matter effects. 




(5.2) 
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FIG. 3: 90%, 95% and 99% CL contours for the setup III. The left column is simulated for ^13 = 1° 
and 5 = 90° assuming normal mass ordering and 623 = 45°. The right column is the same but 
for = 3°. The upper row is the contribution of the beta-beam, the middle row is the electron 
capture channel with the lower row being the combination. 
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FIG. 4: 90%, 95% and 99% CL contours for setup III-WC with solutions from discrete degeneracies 
included for ^13 = 1° (left panel) and ^13 = 3° (right panel) for different values of the CP-phase, 

5 = -90°, 0°, 90°, 180°. 

As for the previous case, the minimum of the continuum solutions for 613 is located at 

6 = — A13L/2 = —tt/2] shown clearly in Fig. [31 

The power of the combination of the beta-beam and electron-capture channels is in the 
difference in phase and in amplitude between the two fake sinusoidal solutions; their combi- 
nation selects a narrow allowed region in the parameter space, much more constrained then 
the two separate techniques. This effect is clearly visible in the lower plots in Fig. [31 

For ^^13 = 3°, there is still some intrinsic degeneracy that cannot be completely removed at 
99% CL. The marked difference between the beta-beam alone and the combination with the 
electron capture in this case demonstrates the importance of data from the high energies. 
For sin^ 26'i3 ~ 10^^ and 6 ~ 7r/2, the beta-beam configuration is able to determine the 
allowed region in the {613-6) parameter space with relatively good accuracy. For this range 
of 613, the dominant interference term helps in resolving any degeneracy. This is not the 
case for other values of 6. 

In Fig. [H we show the results for setup III-WC, with the effects of the hierarchy clone 
solution taken into account. From a comparison of Figs. [3l and [H the increase in event 
rates improves the results substantially. However, owing to the relatively short distance, 
L = 650 km, the mass ordering can be determined only for large values of the mixing angle 
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6'i3 (see below). The hierarchy degeneracy worsens the abihty to measure ^13 and S with 
good precision, especially for negative true values of 6. 

Setups IV and IV- WC 

The use of Boulby instead of Canfranc for a boost of 7 = 369 serves as an intermediate 
case with respect to the last two setups in the sense that the position of the electron capture 
is neither on first oscillation maximum or far into the tail. Consequently, the beta-beam 
has some, but not a considerable coverage of the first oscillation region. Boulby provides 
a much longer baseline, L = 1050 km. This has two contrasting effects on the sensitivity 
to measure CP-violation: on one side it provides sufficient matter effects to resolve the 
hierarchy degeneracy even for small values of ^13; on the other, it decreases the available 
statistics with respect to Canfranc. 

In Fig. [5] and Fig. [6], we report the results for setup IV and setup IV- WC, respectively, 
for normal mass hierarchy and for the cases 613 = 1° (left panels) and 61^ = 3° (right 
panels). In Fig. [5], we show contours for 6 = 90°. Similarly to Figs. [2] and [3l in Fig. [5] 
we present the results for the /?'^-decay flux only (upper row), electron capture flux only 
(middle row) and their combination (lower row). Fig. [6] shows the results for four values 
for 6, including the hierarchy clone solution (blue dashed contours). The electron capture 
channel displays similar behaviour to the Canfranc high boost case (setup III, middle row 
of Fig. [3]) at ^^13 = 3°. The minimum value of the continuum solutions of 61^ is located at 
6 = —66° and the allowed range of values is ~ 2°. For the beta-beam (upper rows), the 
first oscillation maximum energy bin contributes only marginally to the overall sensitivity, 
as can be understood from Fig. [TJ Therefore, the allowed region in Fig. [5] has a different 
shape with respect to the case shown in Fig. [3l In addition, the smaller count rate results in 
a poorer resolution. However, the synergy between beta-beam and electron capture is more 
important here. The lack of concurrence of the beta-beam and electron capture allowed 
regions implies that their combination constrains 6 and ^13 in small ranges. This can be 
clearly seen in Fig. [51 

In addition, the longer baseline allows for a good determination of the mass ordering 
(see below), eliminating more degenerate solutions and providing an improved sensitivity to 
CP-violation with respect to setups III and III-WC. Comparing Fig. [6] with the lower row 
of Fig. [5l the improvement from setup IV to setup IV- WC is noticeable. 
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FIG. 5: Same as Fig. [3] but for setup IV. 
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FIG. 6: Same as Fig. H but for setup IV- WC. 
B. CP-violation discovery potential 

We now consider the CP-discovery potential for the various setups. For the low-7 options, 
setups I and II, the sensitivity is very limited, in agreement with the findings already reported 
in Fig. [21 Henceforth, we will not show results of these two setup as they always possess 
worse physics reach compared to the other setups. 

Setups III, III-WC, IV and IV- WC have a much better physics reach, as shown in Fig. [7] 
and Fig. [HI where the CP- violation discovery potential at 99% CL for the 50 kton TASD or 
LAr detector (setups III and IV) and 0.5 Mton WC detector (setups III-WC and IV- WC) 
are depicted, respectively. In both figures, the CERN-Canfranc baseline is displayed in the 
left panel and the CERN-Boulby baseline in the right panel. In each panel, we present the 
results for the beta-beam only (blue dotted lines) and the combination with the electron 
capture result (red solid lines), both without (thin lines) and with (thick lines) taking the 
hierarchy degeneracy into account. In all cases we note that the addition of the EC channel 
weakens the impact of the intrinsic degeneracy, significantly improving the sensitivity. The 
CP-discovery potential depends on various factors, mainly the available count rate and 
the presence of the hierarchy clone degeneracy. The count rates are important, as can be 
understood by comparing the results for the setups III with III-WC and IV with IV- WC. 
However, it should be pointed out that the WC detector is not optimised for the high 
energies, where the QE cross section is small and the information on the energy for the total 
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FIG. 7: CP-violation discovery potential at 99% CL for setup III (left panel) and IV (right panel). 
In each case, we present the results for the beta-beam only (blue dotted lines) and the combination 
with the electron capture result (red solid lines), both without (thin lines) and with (thick lines) 
taking the hierarchy degeneracy into account. 

cross section is poor. We could expect a similar sensitivity to CP-violation for a TASD or 
LAr detector with exposure a factor only a few times larger than the one considered in this 
analysis. 

The effects of the hierarchy degeneracy are important, significantly more so for 6 < 0, 
where there is a loss of sensitivity to CP-violation by a couple of orders of magnitude in 
sin^26'i3. We can understand this effect by looking at Fig. |H For the shorter CERN to 
Canfranc baseline, we note that for positive 6 either the hierarchy can be determined or, 
where not, the hierarchy clone solution significantly overlaps with the true one. In the case 
of negative values of 6, the hierarchy cannot be resolved even for large values of ^13; the 
clone solution stretches into CP-conserving values preventing the possibility to determine 
CP- violation. The inclusion of the high energy EC channel helps in resolving the degeneracy, 
with CP-discovery down to sin^ 2^13 ~ 3 x 10"^ at 99% CL (left panel of Fig. E]). 

The CERN to Boulby baseline (setups IV and IV- WC) has stronger degenerate effects, 
but it also provides a better ability to resolve them. The results with and without the 
hierarchy clone solution are not significantly different since, for the values for which one 
has sensitivity to CP-violation, the hierarchy can be resolved. For setup IV, the reach is 
limited to large values of sin^ 26'i3. Using the WC detector, the much larger count rate brings 
significantly better results: CP-violation can be established for a large fraction of 6 values. 
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FIG. 8: Same as Fig. [7] but for setup III-WC (left panel) and setup IV- WC (right panel). 

even for sin^ 2^13 ~ few x 10"^ at 99% CL. 

Comparing the two locations of the detector, we notice that the shorter baseline (CERN- 
Canfranc) has a slightly (significantly) better reach for CP-violation at positive (negative) 
values of 6 than the longer baseline (CERN-Boulby). The longer option, however, performs 
slightly better at negative S if the hierarchy is known to be normal and significantly better if 
the ordering is not determined. This is because the longer baseline can identify the neutrino 
mass hierarchy for these values of ^13, therefore resolving this degeneracy. 

C. Mass hierarchy determination 

In Fig. [9l we present the results for the neutrino mass hierarchy determination, but 
only for the setups with a 0.5 Mton WC detector (setups III-WC and IV- WC). We do not 
consider the CERN-Frejus cases; the shorter baseline being unable to distinguish the type 
of hierarchy. 

In both cases, the contribution from the beta-beam channel is shown in blue dashed lines 
and the result for the combination with the electron capture channel is shown by the red solid 
lines. As matter effects are more important at high energies, we see that the inclusion of the 
electron capture flux improves the results, and in particular for the low values of sin^ 2^13 for 
which the measurement is possible. However, the chances to determine the mass hierarchy 
are very limited for the CERN-Canfranc baseline, never reaching more than 30% of the 
values of the CP-violating phase 6. On the other hand, the CERN-Boulby baseline, with its 
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FIG. 9: Fraction of 5 for which the neutrino mass hierarchy can be determined at 99% CL for 
setup III-WC (left panel) and IV- WC (right panel). In each case, we present the results for the 
beta-beam only (blue dotted lines) and the combination with the electron capture result (red solid 
lines) . 

larger matter effect, represents a much more promising setup, for which the determination 
of the mass hierarchy would be possible for all values of 5 for sin^ 2^13 ~ few x 10^^, and 
with a 50% probability for sin^ 29i^ ~ few x 10^^. 

VI. DISCUSSION 

Let us note that by the time this experiment could possibly take place, there will be much 
better knowledge of the neutrino oscillation parameters, improving quite considerably the 
results presented here. Throughout this study, we have adopted a very conservative stance, 
namely the assumption that the errors on the neutrino parameters will remain the same by 
the time this experiment might start taking data. However, this is very likely not to be the 
case. Here, we present the results when no errors are included on the assumed values of 
the neutrino oscillation parameters, i.e., they are known with perfect precision. The actual 
performance of the experiment would lie in between these results and those presented in 
previous sections. For brevity, we only consider results for setups III-WC and IV- WC. 

The CP-discovery potential is shown in Fig. [101 where the left (right) panel represents 
the setup III-WC (setup IV- WC) and the beta-beam only (beta-beam plus electron capture) 
performance is shown by the dashed blue lines (red solid lines) . From Figs. [8] and [TOl the 
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FIG. 10: Same as Fig. [SJ but with negligible error in the value of the assumed neutrino parameters. 

difference between use of the present uncertainty in the neutrino oscillation parameters and 
no uncertainty is seen to be small. 

On the other hand, we show in Fig. [TT] the extent to which mass hierarchy determination 
could improve with better knowledge of the assumed neutrino oscillation parameters. Again, 
we only show the case of setup III-WC (left panel) and setup IV- WC (right panel) with the 
same designations as previous, i.e., blue dashed lines for the beta-beam only contribution 
and red solid lines for the performance of the total flux (i.e., adding the electron capture 
flux). The difference between Figs. |9]and[TT]is substantial. While the qualitative behaviour 
of the relative contribution of the beta-beam only part is very similar, in the case of perfect 
knowledge of the assumed parameters, the determination of the neutrino mass hierarchy is 
possible for 50% of the values of 6 down to sin^ 29i3 8x 10^^ for setup IV- WC (right panel 
of Fig. [TTll . This improvement represents about an order of magnitude with respect to the 
case depicted in the right panel of Fig. [9l 



VII. SUMMARY AND CONCLUSIONS 



Determining the value of ^13, the type of neutrino mass ordering and the presence of 
leptonic CP-violation will be one of the main goals of the future experimental neutrino 
program which is under intensive discussion at present. In the present article we have 
studied a new type of experimental setup which combines a beta-beam with an electron 
capture beam. This can be achieved naturally by using nuclei which can decay into both 
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FIG. 11: Same as Fig. [21 but with negligible error in the value of the assumed neutrino parameters. 

channels. We have studied this idea using the nuclide ^^^Yb which has favourable beta- 
decay and electron capture branching ratios, and only a small alpha decay contribution. 
This combination is very powerful as the EC channel provides a high energy signal at a 
well known energy, while the beta-beam provides coverage of the first and second oscillation 
maxima. The allowed regions in the (^13, 6) plane for the two separate channels have a limited 
overlap resulting in a good resolution of the intrinsic degeneracy. We have understood the 
main features of this synergy by an analytical study of the oscillation probability. It should 
be stressed that this setup does not require two polarities but reaches a very good sensitivity 
by only using the neutrino channel through full exploitation of the oscillatory pattern of the 
appearance probability. 

We have performed a detailed study of the dependence of the physics reach of this ex- 
perimental technique by considering six different setups: two values for the ion boost factor 
7 = 166 and 369; two choices for the detector: a 50 kton TASD or LiAr and a 0.5 Mton 
WC detector; and three baselines: CERN-Frejus, CERN-Canfranc, CERN-Boulby. This 
allowed us to study the impact of the count rate, choice of baseline and the tuning of the 
energy of the beta-beam and EC beam to the oscillatory pattern. We find that the setups 
with low gamma and 50 kton detectors have very poor physics reach, owing to the limited 
event rates. The information on CP-violation is mainly provided by the high energy EC 
signal. We studied the options with 7 = 369, the highest value of the boost factor allowed 
by an upgraded SPS. Setups III and III-WC, which use the CERN-Canfranc baseline, have 
larger count rates and a better tuning of the beam to the oscillatory pattern, with respect 
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to their CERN-Boulby counterparts: setups IV and IV- WC. This results in a very good 
abihty to measure the parameters, see Fig. [HI In particular these setups provide the best 
sensitivity to CP-violation for positive values of 6. However, for negative 6, owing to the 
relatively short distance, the type of hierarchy can be resolved only for very large values of 
^13. The sign- degeneracy prevents discovery of CP- violation in this case, see Fig. HJ The 
CERN-Boulby setups, IV and IV- WC, suffer from smaller count rates and poor tuning of 
the beta-beam to the oscillation pattern. However, they provide a much better determina- 
tion of the hierarchy and possess a good reach to CP-violation for 6 < 0, even if the mass 
ordering is not known. Comparing the two baseline options, if the hierarchy is known to 
be normal from other neutrino experiments, the CERN-Canfranc option has an improved 
physics reach, while if the ordering is not known, the CERN-Boulby baseline outperforms 
the shorter option. For the high statistics scenario, one gets sensitivity to CP-violation 
down to values of sin^ 26'i3 ~ 3 x 10^^ at 99% CL for a WC detector at Canfranc, and 
sin^ 26*13 ~ 10"^ for a WC detector at Boulby. 

In conclusion, we have presented the novel idea of using a single beam which combines 
neutrinos from beta and electron capture decays and have demonstrated the physics reach 
of several possible setups. We have shown that the combination of these two types of beams 
achieves remarkable results. This could naturally be done with the use ^^^Yb, which has 
comparable beta-beam and electron capture branching ratios. As both beams are produced 
from a single isotope, this combination cannot be further optimised. An analogous setup 
would be obtained if a beta-beam and an electron capture beam sourced from different ions 
are combined. In this case, further optimisation of the experiment would be allowed, for 
suitable choices of baselines, Lorentz boost factors, detector size and technology, possibly 
achieving an even better physics reach than the one found in the present study. 
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